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Abstract: We introduce a new paradigm in Composite Dark Sectors, where the full Stan¬ 
dard Model (including the Higgs boson) is extended with a strongly-interacting composite 
sector with global symmetry group Q spontaneously broken to 77 C Q. We show that, 
under well-motivated conditions, the lightest neutral pseudo Nambu-Goldstone bosons are 
natural dark matter candidates for they are protected by a parity symmetry not even bro¬ 
ken in the electroweak phase. These models are characterized by only two free parameters, 
namely the typical coupling go and the scale fn of the composite sector, and are there¬ 
fore very predictive. We consider in detail two minimal scenarios, SU{3)/[SU{2) x t/(l)] 
and [SU{2)'^ x U{1)\/[SU{2) x ?7(1)], which provide a dynamical realization of the Inert 
Doublet and Triplet models, respectively. We show that the radiatively-induced potential 
can be computed in a five-dimensional description with modified boundary conditions with 
respect to Composite Higgs models. Finally, the dark matter candidates are shown to be 
compatible, in a large region of the parameter space, with current bounds from dark matter 
searches as well as electroweak and collider constraints on new resonances. 
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1 Introduction 

The Standard Model (SM) of elementary particles has been tested to an impressive level 
of accuracy in a large amount of experiments. In particular, the current Large Hadron 
Collider (LHC) data do not reveal any significant departure from the SM predictions, 
leaving the experimental evidence of new physics mainly in the non-vanishing neutrino 
masses and the Dark Matter (DM) observations. The concrete nature of the new physics is 
however still a matter of study. Nevertheless, concerning DM, a weak interacting massive 
particle (WIMP) has long been a prime candidate. Furthermore, in light of the celebrated 
experimentally established Higgs sector, scalar WIMPs have earned some attention. The 
addition of a singlet elementary scalar, for instance, has been detailed considered in the 
recent years (see for example Refs. [1-5]). More elaborated models comprise extensions 
with composite scalars. Among them, composite Higgs models (CHMs) deserve special 
attention, inasmuch as they provide an appealing alternative to supersymmetry. In CHMs 
with DM candidates, both the Higgs and the WIMP particles are assumed to be pseudo- 
Nambu-Goldstone bosons (pNGBs) of a new strongly-interacting sector, with a global 
symmetry breaking pattern 0/71. Hence, they can be naturally much lighter than the 
new physics compositeness scale. Explicit models have been worked out for instance in 
Refs. [6-9]. 

In this article we want to address the phenomenological implications of a related but 
different scenario. We consider SM extensions in which only the DM particles are in the 
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pNGB composite sector, being the SM matter content (including the Higgs boson, H) fully 
elementary. (Similar approaches have been previously considered in the literature. See 
e.g. [10-13]). The global symmetry of the composite sector is therefore only broken by 
the SM gauge interactions. As we shall detail below, this class of models presents several 
interesting features: if the coset is symmetric (see footnote 2) the DM candidate is 

protected by a parity symmetry, not even broken by the loop-induced Coleman-Weinberg 
potential; (ii) this class of models turns out to be extremely predictive, for there are only 
two free parameters involved, corresponding to the scale of compositeness fa and the 
typical coupling of the DM sector; and (in) in the simplest scenarios, the values of go 
and /d reproducing the relic abundance make the model evade current constraints while 
being under the reach of future experiments. In the following, we will discuss the general 
structure of these models, while detailed calculations will be provided for two minimal 
cosets, namely [517(2)2 x U{l)]/[SU{2) x [7(1)] and SU{3)/[SU{2) x [7(1)]. 

The article is structured as follows. In section 2 we describe the general properties 
of this scenario, and classify the smallest realizations that can be considered distinguish¬ 
ing those phenomenologically viable. Based on purely symmetry and scale arguments, we 
estimate the size of the couplings, masses and hence the expected phenomenology on the 
minimal setups. We write the complete phenomenologically relevant Lagrangian in sec¬ 
tion 3. In that section the Coleman-Weinberg potential is also computed in a holographic 
framework. Special attention is paid to the way in which this formalism, originally describ¬ 
ing CHMs, is redefined. Indeed, this new setup requires the boundary conditions in the 
ultraviolet (UV) brane to be properly modified. In section 4, we consider current bounds. 
We show that these are dominated by relic-abundance measurements and LHC searches of 
long-lived charged particles and heavy resonances decaying into pairs of SM fermions. We 
conclude in section 5. 

2 A different look at Composite Models 

CHMs [14-18] were originally proposed as a compelling solution to the hierarchy problem. 
The Higgs boson arises as a bound state of a new strongly-interacting sector, with a global 
symmetry group G spontaneously broken to Ti C Q. Thus, its mass is protected by its 
finite size. The Higgs boson is assumed to be a (pNCB) of the global symmetry breaking 
pattern. Hence, the Higgs boson mass can be naturally at the electroweak (EW) scale if the 
new physics scale / is around the TeV. The symmetry G is explicitly broken in two ways: 
by the gauging of only the SM gauge subgroup ^SM of 71, and by the linear mixing of the 
composite resonances with the elementary SM fields. As a consequence, the physical fields 
are admixtures of composite and elementary states (partial compositeness [19]). Then, 
the larger the mixing, the larger the interaction of an elementary field with the strong 
sector —in particular with the Higgs boson, which is fully composite— and hence its mass. 
Thus, contrary to what happens in the SM, the Higgs potential is dynamically generated. 
However, as hrst stated by Witten [20], the radiative contribution from gauge fields generate 
a potential whose vacuum expectation value (VEV) is aligned in the direction that preserves 
the gauge symmetry. Hence, the linear mixing between the elementary fermions and the 
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Figure 1. Pictorial representation of a Composite Dark Sector. The SM gauge bosons are the only 
responsible for the interaction between the elementary sector (including the Higgs boson) and the 
composite sector. 

composite resonances is not only a requirement to correctly reproduce the fermion masses, 
but also to achieve a realistic EWSB pattern. However, both EW precision data (EWPD) 
[21] and the recent LHC searches [22] are getting more and more in conflict with natural 
arguments [23-28]. Indeed, current analyses are already excluding the verge of the natural 
partner parameter space. ^ Our aim here is to point out that, in light of the previous 
discussion, this mechanism can be used not to provide a realistic CHM, but a viable DM 
explanation. To be concrete, we extend the SM with a strongly-interacting sector with 
symmetric ^ coset Q j'H. The SM, including the Higgs boson, is considered to be completely 
elementary, while we assume that only the SM gauge interactions are the responsible for 
the breaking of the global symmetry in the composite sector. We will prove that the 
neutral pNGBs in this scenario (see figure 1) are natural DM candidates. Indeed, if Q jl-L 
is symmetric, only terms with an even number of pNGBs are present in the non-linear 
description of the new strong sector. That means that the composite sector automatically 
respects a Z 2 symmetry for which the pNGBs vr transform as vr —>■ —vr. Glearly, the 
gauging of ^SM ^ T-i Q breaks explicitly the global symmetry, but keeps Z 2 conserved. 
Hence, given that fermions are assumed to be fully elementary, the parity symmetry is 
not explicitly broken. Not even spontaneously broken, inasmuch as the VEV in the loop¬ 
generated Goleman-Weinberg potential is expected to live in the EW-preserving direction, 
according to the Witten result mentioned above. Furthermore, the pNGB masses are not 
restricted by experimental searches to be as low as in the Higgs case, what allows the 
composite resonances to be heavy enough not to conflict with current constraints. As a 
matter of fact, non-singlet scalar DM particles have been found to predict the correct relic 
abundance specially for large masses (see for example Refs. [29-31]). These results will be 
explicitly stated in the next section. All these features make this framework a promising 
scenario for DM. Before going further, however, let us classify the possible models fitting 

^Recently, different models have been proposed to sensibly relax this tension. See for instance Ref. [28]. 

^By symmetric coset, we mean that in which the broken generators, X, commute as [X*, X^] = ifijkT’^, 
where T stands for the unbroken generators. 
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these requirements, from the simplest to more involved realizations. ^ 

• T-L = SU{2) X 17(1), Q = SU{2) x 17(1)^. In this case, the broken 17(1) commutes with 
the whole SM gauge group. Thus, gauge-boson loops do not generate a potential for 
the corresponding neutral pNGB, for it is not charged under the SM gauge symmetry. 
As a consequence, it is an exact NGB and hence massless. This simplest model should 
be therefore disregarded from the phenomenological point of view. 

• % = SU{2) X 17(1), G = SU{2) X 17(1)"'. For any n, the argument above still applies. 

• T-L = SU{2) X 17(1), G = 517(2)^. In this case, only two charged states emerge. In 
light of the strong bounds on stable charged particles [31-34], this kind of models 
are not phenomenologically viable. More in general, pNGB spectra with only neutral 
fields are not expected. Consequently, the minimal realistic models do have at least 
three degrees of freedom, corresponding to one neutral and two charged states. 

• T-L = 517(2) X 17(1), G = 517(2)^ x 17(1). This is the simplest realistic realization. A 
real scalar triplet with hypercharge T = 0 appears in the spectrum. We will analyze 
this model in detail in next sections. 

• Ti = 517(2) X 17(1), G = 517(2)^ x 17(1)". As in the first case, massless scalars are 
expected. Thus, we no longer consider this coset. Besides, further extensions with 
more 517(2) insertions give rise to coset spaces larger than the ones considered below. 

• 7i = 517(2) X 17(1), G = 517(3). This symmetry-breaking pattern has been previously 
considered in the CHM literature (see for instance Ref. [35]). However, it has been 
disregarded in the recent years for it gives large corrections to the p ~ 1 parameter 
(it does break the custodial symmetry in the Higgs sector). In our setup, instead, 
the Higgs boson lives in the elementary sector and thus we will also take this coset 
into account. 

Larger group combinations either provide cosets with more degrees of freedom or scenarios 
phenomenologically similar to the ones we have enumerated previously. A prime example is 
50(5)/50(4). This coset structure stands for the minimal GHM and has been largely stud¬ 
ied in the literature [36]. This model has the advantage of containing no anomalous repre¬ 
sentations, that potentially distort the parity symmetry [7]. At any rate, in the following we 
only consider in detail the cosets 5f7(3)/[5f7(2) x [/(I)] and [517(2)^ x 17(1)]/[517(2) x 17(1)], 
assuming that their UV completions make them anomaly-free. Finally, a last comment is 
in order. Dark multiplets with fractional hypercharge have been typically ignored so far. 
The reason is that these multiplets are always complex, and so two neutral states are pre¬ 
dicted. The problem arises because these states are generically degenerated in mass, and 
they are hence strongly constrained by direct detection experiments [37]. The exception is 

^Hereafter, both H and Q are always assumed to be multiplied by the unbroken color group SU (3) 
too, which is only omitted in the equations for a simpler reading. Although its inclusion could be also 
disregarded, we decided to take it into account since it just leads to more conservative bounds from direct 
detection. This argument also applies to other spectator groups to first approximation. 
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the well-known case where a second EW doublet <I> is added, the so-called Inert Doublet 
Model (IDM) [38-40]. In this case, the operator A [(fl^<l>)^-|-h.c.] can be written at the tree 
level. After EWSB, this operator introduces a small splitting in the mass squared of order 
~ between the two neutral components. This operator, however, does not arise at the 
quantum level in the Composite Dark Sector scenario provided by SU{3)/[SU{2) x 1/(1)]. 
The reason is that, in the massless limit, the pNGB sector respects a 17(1) symmetry (that 
includes the Z 2 parity) under which —)• exp (i0)4>, being this symmetry still exact after 
gauging the SM gauge group. At any rate, it can be always assumed this symmetry to 
be broken at a higher scale. Not only in the Y = 1/2 case, but also for larger quantum 
numbers. As a matter of fact, for any SU{2)l multiplet cj) with half-integer hypercharge 
and isospin I satisfying the relation Y = I, this splitting can be achieved by a higher¬ 
dimensional operator with m = 21 Higgs insertions and an arbitrary number n of Higgs 
singlet operators X[{H^H)^{HH — ■ H+ h.c.], where H = ia^H*. Note that a new 
physics scale is anyway required in order to stabilize the Higgs mass. Hereafter we assume 
that this is the case in the coset SU{3)/[SU{2) x 17(1)]. Provided the splitting is small 
(as it can be expected if the additional 17(1) symmetry is broken at a large scale), the 
dependence of other constraints as well as future searches on it is negligible. Thus, the 
model is effectively still described by just two parameters. 

A first good approximation to the size of the masses and couplings of the degrees of 
freedom that arise in both [517(2)^ x 1/(1)]/[517(2) x 17(1)] and 517(3)/[517(2) x 17(1)] can be 
obtained considering the symmetries and scales involved in this framework. Indeed, given 
that the new sector is assumed to be strongly coupled, go is expected to be 1 < < dvr. 

Thus, the new resonances are expected to have a mass rrip ~ gofo ~ few TeV in order not 
to conflict with current constraints. Thus, the pNGBs vr (which are massive only because 
of the explicit breaking of the global symmetry) are naturally expected to live around (or 
slightly below) ~ 1 TeV. In addition, the non-derivative scalar interactions, being generated 
by loops of SM gauge bosons, are predicted to be sub-dominant with respect to the gauge 
interactions. Besides, as both the electromagnetic charged and neutral pNGBs come in 
complete representations of the SM gauge group, the mass difference among themselves 
can only arise in the EW broken phase, being hence proportional to some power of the 
EW VEV V ~ 246 GeV. (As a matter of fact, there is only one renormalizable operator 
that can break the degeneracy after EWSB, given by (7llr“i7)Mab(7rlr^7r), where r and T 
are the Pauli matrices and the 517(2) generators in the representation vr belongs to in the 
spherical basis, ^ and Mab is the three times three anti-diagonal matrix adiag(l, —1,1).) All 
together implies a small splitting between the neutral and the charged states, that makes 
the latter typically long-lived. As we shall discuss in section 4, searches of these particles 
in the current LHG data [41, 42] provide one of the strongest constraints in some scenarios. 

^The Pauli matrices in the spherical basis are written as = ±(cr^ =pi(T^)/2 and Y = cY j"I. The same 
applies for larger representations. 


- 5 - 



3 Gauge interactions, scalar potential and masses 


Let us consider the two cosets [SU{2)‘^ x U{l)]/[SU{2) x t/(l)] and SU{3)/[SU{2) x ?7(1)]. 
The group matrices are explicitly written in the Appendices A and B respectively. The 
dynamics of the composite DM boson at low energy is described by a non-linear sigma 
model, with non-linear interactions parameterized by the scale Jd- Following the CCWZ 
formalism [43, 44], the derivative and gauge interactions among the pNGBS are given by 
the trace of d^d^, where stands for the perpendicular projection of the Lie-algebra- 
valued Cartan one-form = —iU^D^U = -|- E^T°‘, with U = exp {—i ^ fo), 

X°‘ the broken generators, r“ the unbroken ones and the covariant derivative: 


+ igs^gl + igT^wl + ig'YB^ 


which can be also written as 


Df,ip = 


dfi + igs—g'^ + + T ) H —— (Ts — s'^Q) Zy, -|- ieQA^ 


(3.1) 


(3.2) 


if we define ±iT^ and = {W^ ^iW‘j^l\f2. We have denoted by sw (cry) the 

sine (cosine) of the SM Weinberg angle 9\y and e is the electromagnetic unit charge. For 
[SU{2)‘^ X U{1)]/[SU{2) X f7(l)] the relevant Lagrangian ^ reads: 


C = ff2(vr0)2py+vy/^- + 


igW^^{'K^^T: )— \g^W'^W^^TT vr -|-h.c. 


(3.3) 


+ ff'lF+W^-Tr+Tr- + ^{siy - if ^/^(^r+^vr-) 

Cw 

+ e^A^A^TT~^TT~ -t- ieA^(7r+^7r“) -h — (s^ - l)A^Z'^7r+7r“ 

cw 

+ egA^7r°lF^+7r“-h — (s^ - l)W+Z^7r°7r“-Fh.c. + ^^[d^{Tr^f]d^{TT~^Tr~). 
Cw i ^Jd 


For SU{3)/[SU{2) x [/(!)] it reads 


c= + 

4 ^ 

^2 


+ - ^g^WfW^^+TT-TT- +ll.C. 


(3.4) 


+ -b ^(1 - 2s^)2z^z^7r+7r + n ) 

-b e^A^A^vr+vr” -b ieA/,(7r+^7r“) + -{I - 2s^)A^Z^7r+7r“ 

Cw 

— 4 ttOw+^- _ — w+z^nV" + —+ h.( 


+ ( ^A.n^w+TT- - 
.,2 


2 - M ^ ” A* '■ 9^,,, " 


+ 


g 


2cw 
1 


2cry 


8^2 

ocw 


+ [a^|nO|2]dA‘(7r+7r-), 


2 / 


D 


®For the computation of the potential (see next section) all the terms with larger powers of l/fo are 
included. So, once the whole potential is known, the study of the DM and the LHC phenomenology can be 
completed with the addition of the Lagrangian above. 
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In the equations above, + iA9, and ^ stand for the neutral pNGBs and 

7r+ and 7r“ are the charged ones. In addition, we have defined = fd^g — {d^f)g. 

We have omitted higher-dimensional operators, which are suppressed by larger powers of 
l/fn, with the exception of the last term in eqs. (3.3) and (3.4). The reason is that, 
although a detailed computation of gamma-ray observables is beyond the scope of our 
study, ® it is worth mentioning that these operators can play an important role on such 
processes. Indeed, processes like —)• 'y'y{Z) receive contributions from loop triangles 

with virtual vr^. The coupling between the neutral and the charged states is usually given 
by the quartic term (7r^)^7r'’'7r“ in the potential. However, the contribution from the 
derivative interaction, « can easily exceed the coupling A.n-o^±. Specially in these 

composite models, in which the latter is radiatively induced. 

In order to estimate the strongly-coupled effects, we work out a holographic scenario 
described by a modified five-dimensional (5D) description of CHMs, where the full SM lives 
on the UV brane and just gravity and gauge interactions propagate through the bulk of 
the extra dimension. In particular, we consider an AdSs space with metric [47] 


ds^ = a‘^{z){r]f^ydx^dx'' — dz‘^) 


f— {g^Ax^dx" - dz^), 


(3.5) 


where 2 ; G [R, R'] is the coordinate of the additional spatial dimension and R and R' are 
the positions of the UV and the infra-red (IR) brane, respectively. In order to effectively 
describe the two breaking patterns introduced before, we will extend in the bulk of the extra 
dimension the SU{2)l x U{ 1 )y gauge symmetry of the UV brane to SU{2)i x SU{2)2 x 
U{ 1 )y and SU{3) x U{l)x, respectively, where the additional U{l)x gauge group has been 
added to correctly reproduce the observed Weinberg angle. On the other hand, as the 
IR boundary conditions parametrize the spontaneous symmetry breaking of the Goldstone 
symmetry by the strong dynamics [48, 49], the 5D bulk gauge symmetry will reduce on the 
IR brane to SU{2)i x U(l)y and SU{2)i x U(l)8 x U{l)x, respectively. More specifically, 
we define 

= + = M = g,5, i = l,2,3, (3.6) 

in the [SU{2)‘^ x U{1)]/[SU{2) x 17(1)] case, with and R^j^ being the 5D gauge fields 
associated to SU{2)i and 517(2)2, respectively. For the 517(3)/[517(2) xl7(l)] case, instead, 
we have 


Bm 

C-(j) 


S(j>V3W^ C(j)XM, 
95 


\l9l+9\ 


9x 

— 


\l9l + g\ 


M = /i, 5 , 

(3.7) 


where g^ and gx are the dimensionful 5D gauge couplings of SU{3) and 17(l)x, while 
VF^, I = 1,...,8, and Xm are the corresponding 5D gauge fields. In each case, the 

®As well as other relevant gamma-ray spectral features, see e.g. [45, 46] for the IDM case. 
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(3.8) 


corresponding boundary conditions on the UV before EWSB read ^ 
ie;(+,+), B^{+,+), i = 1,2,3, 

and 


^/i(+)+)) ,+), z —1,2,3, a —4,5,6,7, (3.9) 

respectively. In the previous equations, the first (second) +/— denotes Neumann/Dirichlet 
boundary conditions at the UV (IR) brane and we have omitted those of the four-dimensional 
scalar counterparts (/r —)> 5), which have opposite boundary conditions. For both cosets, 
the gauge fields with UV Neumann boundary conditions, RV and will be associated 
with the SM-like EW gauge bosons, whereas the scalar components of the ones having 
Dirichlet boundary conditions at both branes will provide the corresponding Goldstone 
degrees of freedom 

Xiix,z) = f{z)7r\x) + ..., TU“+3 =/(z)7 r“(x) + ..., i = l,2,3, a = 1,2, 3,4, (3.10) 


where the dots stand for non-physical Kaluza-Klein (KK) resonances and [35] 


f{z) = a (z) 


rR' 

/ dz' a-\z') 

Jr 


- 1/2 


(3.11) 


Since in these scenarios, and contrary to the usual CHMs, the Higgs doublet is localized 
on the UV brane, the boundary conditions after EWSB at the UV brane for the gauge 
bosons with Neumann boundary conditions will be modified. It is thus convenient to 
define the usual physical combinations = (lU^ =F ^H^m)/\/ 2 and 


Am 

cw 


SW^M + CwBm, 
95 

Voi + gf 


Zm 

sw 


CW^M ~ SwBm, 
95 

Vgi + gf 


(3.12) 


with (/g being the dimensionful 5D gauge coupling associated to the hypercharge U{1)y, 
given by 


95 


95gx 

+ 9x 


(3.13) 


in the SU{3) case. ® Therefore, the UV boundary conditions for and after EWSB 
read 



(3.14) 

(3.15) 


^We will assume for the moment that the EWSB pattern is the usual one, checking explicitly later that 
this is indeed the case. 

®For the holographic description of SU{2)i x 517(2)2 —>■ SU{2)l we have assumed a Plr symmetry 
SU{2)i 517 ( 2)2 within the composite sector, leading thus to = g^^^ = V2gs and reducing by one the 
number of 5D input parameters. In this case, g'g is just the 5D gauge coupling associated to 17(l)v. 
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whereas all the rest remain the same. After imposing the UV boundary conditions, this 
leads in particular to the following KK decomposition ® 

\{x,z) = ^a^C{rnl,z)A^if>{x), (3.16) 


Z^{x,z) 

Z'^{x,z) 


n 

G{m^,z) + -^^{gl+g'i)S{m 
rrin ^ 

^ z) 

n'> ^ ) 

Zl”>(rc) 


G{m^ ,z)+ g1S{m'Y , z) 

irin ^ 

lU±W(x), 





n 


(3.17) 

(3.18) 

(3.19) 


where C{m, z) and S{m, z) are functions satisfying the bulk equations of motion 


[a{z)m‘^ + dza{z)dz] f{z) = 0, 


(3.20) 


and fulfilling boundary conditions C'(m, i?) = 1, dzC{m,R) = 0, S{m,R) = 0, dzS{m,R) = 
m. They are given by [50, 51] 

C{m,z) = ^mz[YQ{mR)Ji{mz) — jQ{mR)Yi{mz)] , (3-21) 

S{m,z) = ^mz[Yi{mz)Ji{mR) — Ji{mz)Yi{mR)] . (3.22) 

It is possible to express three of the five input parameters in these 5D holographic 
descriptions {v, g^, R, R'} as a function of the other two by matching Gp, mz and 
aem("iz) = e(m^)^/47r to their best SM fit values. The first of these conditions just yields. 


where Gj^\z, z') is the 5D W propagator evaluated at zero momentum, i.e.. 


(3.23) 


Gr(.,.o=-E 




n=l 


m. 


W2 


glv^ 


min( 2 ;, z'Y — R 


(3.24) 


Naively, as usual in these scenarios, one might have expected that the contribution of the 
KK tower would require a shift of the Higgs VEV to absorb it, i.e., v Y However, as 
one can readily see using the holographic basis, all the effects arise from the UV localized 
operator {D^HY, leading to 


1 _ Gf _ 

2u2j^ V2 2u2 ’ 

®Obviously, eq. (3.19) just holds for the SU{3) case. 


(3.25) 
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and thus to u « 246 GeV. It is possible to get analytical expressions for the other two 
constraints by using the approximate expressions of C{m, z) and S{m, z) for mz < mR' <C 1 


C{m,z) K, 1 — / dzia ^( 2 : 1 ) / d 2;2 0 ( 22 ) 


'R 


'R 




rz2 




+ m I dzia {zi) / dz 2 0 ( 22 ) / dz^a ( 23 ) J dz^ 0 ( 24 ) 


JR 

= 1 — -m^R^ 


IR 




+ 


R 

l + 4(^) -5(4) +4(^1 ( 2 +(^) )log(^ 


(3.26) 


R 


RJ 


R 


R 


R 


rz i-z i-zi |^Z 2 

S{m,z) Ri m / d 2 ;ia“^( 2 ;i) — / d 2 ;io“^(zi) / d 2;2 o-iz 2 ) / dz 3 o“^( 2 ; 3 ) 

Jr Jr Jr Jr 


IR 
9 

mz 


2R IGR 




which leads to 


1 55 

mz ^ 


2 7J/2 


1 - 


1 gi v^R 


2 y/RL [ 32 c^i?L L 

where we have defined for convenience the volume factor 

L = log(ii7i?). 

Therefore, one can write at leading order in v‘^R'^ 


and e = 


g5Sw 

Vrl' 


9D 


= g^R-R'^ ^gVL[l + 


1 mL i?'2 1 




2 "2 7 1 ’ ~ sin 20w ( 1 — 

' — S^r L 


1 

8 L 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


where go is the dimensionless coupling in the strong sector. Moreover, since we are no 
longer trying to address the gauge hierarchy problem, there is no need to require 1/7? to 
be roughly given by the Planck scale 1/7? ~ Mpi. Thus, we can consider both 7? and R! as 
free input parameters, or equivalently, go and the scale of compositeness 


Id 


20(1) 

55 


rR' 

/ dz a ^{z) 

Jr 


- 1/2 


2^WV2 

goR' 


(3.31) 


In the above equation and hereafter 2° correspond to the SU{2)‘^ x 7/(1)]/[57/(2) x 7/(1)] 
coset whereas 2 ^ will refer to the 57/(3)/[57/(2) x 7/(1)] one. Since fo is not a pure physical 
quantity and its definition is always arbitrary, we have chosen fo as the scale appearing 
in the argument of the tower contribution to the corresponding pNGB potentials (see 
below). This means in particular that the mass scale of the composite vector resonances 
rup will scale differently with fogo for each coset (see e.g. figure 3). 

In the UV unitary gauge, where all the NGBs of the SM are gauged away, a scalar 
potential l/(/i, tt**^), with k = i,a, will be generated at the quantum level through the 
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interaction of the NGBs with the towers of resonances associated to the 5D gauge bosons 
A^(x, z) = Af,(x, z), Zf,{x, z), Wij^{x, z) [52, 53], 


V(h,7r^} 



n=l 


d'^p 

(W 


log 


p^+ml(h,7r'^ 


(3.32) 


where vr®*'), n G N, are the masses of all possible KK resonances A^^\x) in the 

presence of the background fields vr*^ and h. The previous infinite sum can be exchanged 
by an integral on the Minkowski space [50, 54], 

3 

d/c A;3logp;,^^k(-A:2), (3.33) 

where k = and w G C, is some spectral function, holomorphic in the 

Re(a;) > 0 part of the complex plane and with roots in the real axis encoding the spectrum 
of A^(x, z) in the presence of the background fields and h, i.e., 

Ph,-K'^{'>^n{h,7r^)) = 0, n G N. (3.34) 


Such a function will be proportional to the determinant of the linear system of equations 
resulting from imposing the IR boundary conditions after we remove the NGBs vr*®’ from 
the bulk via the following 5D gauge transformation 

Anix, z) Q.{z)Am{x, z)Q.{z)'^ - (i/g^){dM^{z))9.{z)'^ , (3.35) 


where Am £ Appendices A and B for the definition of the 

different generators) and 

Q{z) = exp ^—ig^X^TT^ J dz'f{z')^. (3.36) 

The terms in the scalar potential V{h,TT^) involving are expected to be finite, due to 
non-locality in the 5D theory [53, 55], as it is manifest by the dependence of the spectral 
function .„.k on the Wilson line 


W = = exp 



exp , 


(3.37) 


which is clearly a non-local object depending on the conformal distance between the branes. 
Even though UV localized terms can give infinite contributions to the V{h) C V{h,7T^) 
potential, they can be renormalized in the usual way. Moreover, any V(h) piece of V{h, tt**^) 
can be shifted to RsM(h), the usual SM potential that needs to be added to (3.33). 

In order to investigate the shape of the pNGB potential V{h,TT^) and to be sure that 
the desired pattern of EWSB is not changed, we perform an expansion of (3.33) in powers 

appears via UV boundary conditions and it leads to analogous expressions to (3.17) and (3.18) with 
the replacement u —> h, since H = + *‘^°] )^- 
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Figure 2. Left) F(/i, 7r°) + Vsm(^) in the [S'C/(2)^ x U{V)\/[SU{2) x t/(l)] coset for = 1 TeV 
and qd = 3.5, where tt^ = 0. Right) V{h,'K°) + Vsm(^) in the SU{3)/[SU{2) x C/(l)] coset for 
/d = 1 TeV and go = 3.5, where = 0. 


of h/fo and sm(n//£)), where IT = obtaining for the [5C/(2)^ x U{1)]/[SU{2) x 

[/(!)] coset 


v{hy) 


Ao + A2 



1 + 


1 

2 


tan^ 6w 


n2 j 



(3.38) 


where vr^ = (vr^ =F and we identify for this coset the neutral state tt^ = We 

have also defined 


3 2yR' 

■ ^ A RC'iVt, R')S{Vt, R')' 

_ _A r dt 

327r2A R'C>{Vi,R'fSiVi,R')' 

3 RS'{Vi,R'f 

3^ A 2R'^ytc'{Vi, R'fs{y R') ■ 


(3.39) 

(3.40) 


(3.41) 


In the above expressions, 5(m, z) and C{m, z) are the Wick-rotated versions of S(m, z) 
and C{m, z), 

C{m, z) = C{im, z) = mz [Ii{mz)KQ{mR) + lQ{mR)Ki{mz)] , (3.42) 

S{m, z) = z) = mz [Ii{mz)Ki{mR) — Ii{mR)Ki{mz)] , (3.43) 

and / denotes the partial derivative with respect to A is an IR cut-off which has 
been introduced to regulate the spurious IR divergences arising from the expansion of 
the Coleman-Weinberg potential [24, 51, 56]. It has been fixed by asking (3.38) to repro¬ 
duce the exact Goldstone mass splitting Am = m^± — m^o. On the other hand, for the 
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SU{3)/[SU{2) X C/(l)] coset, we obtain 


V{h,7r° 


Aq — (7-1-2 sec 






1-1-3 tan' 


9w'^ Ao 


-|- (38 — 20 sec^ 6 w + 12 sec' 


^ 0w^ ^2 


-|- 2 tan 6 w ^2 ( -r~ 


A 

fn 

2 


sin^ ( 2^^ 


fnj 


h y + {A^fY - . 


Id) 


n4 


n 


sin^ ( 2^ h (3-44) 


where now vr^ = (vr^ =p and we have identified for this coset vr^ = vr^ and = tt^. 

The couplings Aq and A 2 are exactly those given by eqs. (3.39) and (3.40) but with a 
different IR cut-off, reproducing the corresponding splitting m^± — m^o = m^± — m^o. 

In order to study the behavior of the scalar potential, we show in hgure 2 the potential 
V{h, 7 r^) + Rsm(^) in both cosets, where is the lightest neutral pNGB and all other 
degrees of freedom have been set to zero. In both cases we have chosen benchmark values 
f£) = l TeV and qd = 3.5. One can readily see from these plots that the interaction with 
the pNGBs does not spoil the Higgs EWSB and that they acquire no VEV as expected. 

We can avoid the use of the IR cut-off introduced before by Taylor expanding (3.33) 
around h = v and = 0 to the renormalizable level. Even though some of the quartic 
self couplings arising in this expansion are still IR sensitive, all the relevant couplings for 
our phenomenological study are IR safe. This is indeed what we have done to obtain all 
relevant couplings henceforth. 

The masses of the neutral and charged states as well as the mass difference are shown 
in figure 3 left and right respectively. The masses of the KK excitations are also shown for 
comparison. As expected, there is a large gap between the scalar and the vector resonances 
due to the pNGB nature of the former. In addition, we observe a rather small splitting 
between the neutral and the charged scalar states. This result is a consequence of the 
potential being loop-induced. Besides, it is worth to point out that the splitting is much 
smaller in [51/(2)^ x U{1)]/[SU{2) x U{1)] than in SU{3)/[SU{2) x H(l)], because it only 
arises at order rn^/ffj in the triplet case. 


4 Current constraints 

The phenomenology of these models is completely described by only two free parameters, 
that we have chosen to be go and fo- We consider a region in the fn — 9 d plane, pa¬ 
rameterized by go £ [1-5,4] and /d £ [1,10] TeV. This is well motivated since, as we will 
see, it contains a large region of parameter space for which the relic density can be exactly 
reproduced. In this region, restrictions coming from EW constraints on the S, T and U 
parameters are not sensitive enough to set important limits on these models, but modifica¬ 
tions of W and Y do impose non-negligible bounds. These are discussed in section 4.1. On 

^^There is a deep relation between this result and the fact that, in the renormalizable triplet scalar model, 
the corresponding quartic coupling does not renormalize under gauge interactions [29]. 
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Figure 3. Top left) Masses of both the charged pNGBs (dashed lines) and the heavier resonances 
(solid lines) as a function of fo for three different values of go in the coset [S'f7(2)^ x C/(l)]/[S'f7(2) x 
f7(l)]. Top right) Same as top left) but for the coset SU{3)/[SU(2) x C/(l)]. Bottom left) Mass 
difference between the neutral and the charged pNGBs as a function of fjj for three diferent values 
of go in the coset [SU{2Y x U{1)\/[SU{2) x f7(l)]. Bottom right) Same as bottom left) but for the 
coset SU{i)/[SU{2) x ;7(1)]. 

another front, Higgs searches are not sensitive to this region of the parameter space, for 
both the pNGBs and the composite resonances are not light enough (see figure 3) neither 
to allow the Higgs decay into invisible particles nor to modify appreciably the Higgs decay 
into SM fields. 

Besides, direct detection experiments and monojet searches provide very weak bounds 
on our parameter space. Indeed, in the first case, nucleon-DM scattering processes are 
mediated either by loop-suppressed processes or by the t-channel exchange of a Higgs 
boson with a strength proportional to the coupling in the quartic //^(tt^)^ term. It has 
been shown in previous Refs. [6] that direct detection experiments are not sensitive to small 
values of this coupling, specially for large DM masses. In particular, values below <0.1 

^^The mediation of a t-channel gauge boson results in a heavier final state, and hence suppressed by the 
small DM velocity. 
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are out of the reach of direct detection experiments for any DM mass. The corresponding 
coupling in our models turns out to be much smaller (~ 10“^) than this value, as can be 
derived from eqs. (3.38) and (3.44). Concerning monojet searches, the rather large pNGB 
masses together with the volume-suppressed couplings of the heavier vector resonances to 
the SM quarks put these models beyond the reach of the current LHC analyses. We have 
explicitly checked this observation implementing the model in MadGraph v5 [57] by means 
of FeynRules v2 [58]. We have generated Monte Carlo monojet events in the parameter 
space region under study. The largest cross sections correspond to the smallest values oi go 
and fo in the SU{3)/[SU{2) x U{1)] model, being of order ss 1 pb. We have subsequently 
passed these events through Pythia v6 [59] and implemented the CMS analysis in Ref. [60] 
(based on an integrated luminosity of 19.7 fb~^ at ^/s = 8 TeV) in MadAnalysis v5 [61] 
for a cut on the missing energy of > 450 GeV. The resulting acceptance times the 
production cross section after the cuts is always smaller than the upper bound on this 
quantity as stated by CMS, namely 7.8 fb at 95% C.L. 

Thus, beyond EWPD, the main constraints come essentially from three sides. First, 
from the measurement of the DM relic abundance Dobs^^ = 0.1199 ±0.0027 [62], which sets 
an upper bound on the contribution for any DM candidate. (These bounds are discussed 
in section 4.2.) Second, from collider searches of long-lived charged particles. Indeed, given 
the small splitting between the neutral and the charged scalars, the latter can be long-lived 
enough to scape detectors, leaving a characteristic trace because of its large mass. We con¬ 
sider these searches in section 4.3. (Note that these particles can be considered long-lived 
for collider experiment purposes, but not at cosmological scales. Thus, cosmological con¬ 
straints as those coming from nucleosynthesis can be neglected [63].) On top of that, there 
is a third source of constraints. These are collider bounds on the heavier resonances, the 
main ones coming from the LHC ATLAS and CMS Collaborations. These last constraints 
are detailed in section 4.4. 


4.1 Constraints from electroweak precision data 

To start with, the T parameter does not receive tree-level corrections from the extra scalars, 
because the EW VEV is aligned in the Higgs direction, as we showed in the previous section. 
Besides, the loop-induced contributions from these scalars to the S, T and U parameters 
are suppressed by the splitting between the charged and the neutral components [64, 65]. 
Hence, given that this is at least at the percent level (see figure 3), these corrections are 
expected to be negligible, in light of the latest measurements: S = —0.03 ± 0.1, T = 
0.01 ± 0.12 and U = 0.05 ± 0.10 [66]. In addition, tree-level corrections to the T and S 
parameters from the heavier resonances are absent. Indeed, this can be easily understood 
in the dual 5D model. If one computes the contributions to the oblique parameters arising 
from the integration of the KK gauge resonances, see e.g. [67, 68], it can be readily seen 
that 


T = 


tan^ 

ALCXprin 


d — 2/3 ± 7 


= 0 , 


s = 


2gjj sin^ 6wv 
LcXori 


2 r 


-/3±7 =0, (4.1) 
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Figure 4. Left) Excluded region in the coset [SU(2)'^ x U{1)]/[SU(2) x 17(1)] in the fn — go 
parameter space. The pink region is excluded by relic density measurements. The largest green 
region is excluded by searches of long-lived charged particles at the LHC. The light shaded region 
is excluded by dijet searches whereas the medium and dark shaded regions are excluded by ti and 
dilepton searches, respectively. The small dashed line encloses the region excluded by EWPD. 
Right) Same as Left) but for the coset SU{3)/[SU{2) x C/(l)]. 


since all the coefficients involved, d, /? and 7 , 

a = RL dz dz'a^{z)[fh{z)fGoiz,z')[fhiz')fa^{z'), (4.2) 

Jr 

/3 = RL dz'GoiR,z')[Uz')fa^iz'), (4.3) 

Jr 

^ = RLGo{R,R), (4.4) 


become equal for a UV localized Higgs {z)[fh{z)]‘^ = d(z — R), where Go{z, z') is the 5D 
Dirichlet propagator before EWSB at zero-momentum after subtracting the corresponding 
zero-mode, 


Gr.{z,z') = 


(1 + 2 log (^)) + 4 (1 + 2 log (f) ) - [i?'2 _ r2 ] ^-1 

ARL 


(4.5) 


and 2 :< = min(z, z'), Zy = max{z, z'). However, we are still left with a contribution to 


^2 ^.2 „.2 

W = Y = 


AL 


AL 


1 1 
Yl 


9D 


fnj 8 ’ 


(4.6) 


arising from four-fermion interactions and which can be in principle non-negligible for 
g'D ~ 1 and not too large values of fo- In order to study the impact of such operators on 
EWPD, we have performed a complete and up-to-date fit to EWPD for 5 = T = 0 and 
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W = Y The allowed values at 95% C.L. are given by the region above the dashed lines 
in figure 4. 

4.2 Constraints from the relic abundance observation 

Assuming the well-motivated standard thermal history for DM, the relic abundance can be 
computed using MicrOMEGAS v4 [71]. Thus, we have also implemented the model interac¬ 
tions in CalcHEP v3 [72] by means of FeynRules v2. We require the computed relic density 
to be < 0.12. The corresponding excluded region in the fn — dD plane is shown in pink 
in figure 4. In the left panel we show the results for the [SU{2)‘^ x U{l)]/[SU{2) x C/(l)] 
case, while in the right panel we consider 517(3)/[517(2) x 17(1)]. Hence, the points in 
the frontier with the white region correspond to the parameter space region in which the 
observed relic density can be explained by a single composite DM particle. Points aside 
this region then require extra degrees of freedom to account for the observed relic density. 


4.3 LHC constraints on long-lived charged particles 

The small splitting between the neutral and the charged states (shown in figure 3) makes 
vr^ long-lived. Indeed, it mainly decays to the vr^ through the emission of an off-shell 
gauge boson, being the total width given by the approximate expression [73] 


g'^a Am® 
dSvr^ m^y ’ 


(4.7) 


where a = 1 (1/2) in the triplet (doublet) case. In [517(2)^ x 17(1)]/[517(2) x 17(1)], the 
lifetime r = l/T is large enough to allow to scape the detectors, if it is produced in high- 
energy collisions. In 517(3)/[517(2) x 11(1)], on the contrary, the decay always takes place 
in the detectors. The decay products are however too soft to trigger the corresponding 
final state in order to make it emerge from the huge W+ jets background. However, in the 
triplet case, the trace of these particles can be still observed for they give rise to anomalous 
energy loss. In fact, there are dedicated analyses to search for this kind of signatures. 
In particular, the CMS Collaboration has reported the latest constraints in Ref. [41]. To 
our knowledge, these are the strongest limits from collider experiments. In that article, 
bounds are provided for different type of charged states. Among them, we can find limits 
for pair-produced staus, which can be directly translated to our model if the theoretical 
cross section for the pair-produced charged scalars is computed. In order this cross section 
to be calculated, we again use MadGraph v5. In figure 4, the green region below the solid 
black line represents the parameter space points for which the cross section exceeds the 
values provided by CMS. 


4.4 LHC constraints on new heavy resonances 

Heavy gauge boson partners are common predictions in models with composite sectors (or 
in their extra-dimensional dual models). In our case, as a consequence of having elementary 

are grateful to Jorge de Bias for providing us the for the EW fit. This fit includes all the 
observables considered in the analysis of [69, 70], updated with the current experimental values. 

^^In the doublet case, could decay also into A° provided that m^o < m^±. 
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SM fermions, these resonances will couple to them with universal and volume-suppressed 
couplings, i.e., ~ g/go- Moreover, for the same reason, the interactions with the SM Higgs 
will be also volume suppressed whereas those to the pNBGs will be on the contrary volume 
enhanced. In particular, this means that unless g/go ~ 1 they will have an important 
fraction of invisible decays making their collider observation extremely challenging even 
for small values of fi). Thus, LHC constraints on these resonances will just be relevant for 
small values of both and fo- Indeed, these vector resonances are expected to be heavy 
and narrow, what allows us to directly translate the limits on new resonances provided 
in the LHC analyses to our model once the corresponding cross sections are computed. 
Among these analyses, we consider searches of: (i) heavy resonances decaying into pair 
of jets. We take the limits from the CMS analysis in Ref. [74], in which the data set 
corresponds to a collected luminosity of 19.7 fb“^ at ^/s = 8 TeV. The points excluded by 
this analysis are given by the light shaded region in figure 4; (ii) searches of tt resonances. 
We focus on the analyses provided by Ref. [75] . The center of mass energy and integrated 
luminosity are the same as before. Civen the small branching ratio into ti, the bounds are 
much weaker. These are shown in figure 4 in the medium dark shaded region; (in) searches 
of dilepton resonances. We use the latest limits provided in Ref. [76]. The corresponding 
excluded region in figure 4 is given by the dark shaded region. There is a last analysis that 
could be potentially useful, namely the search of pair-produced heavy resonances, as those 
provided in [77] and [78]. However, their sensitivity to large resonance masses is still very 
limited, and so the bounds turn out to be completely negligible in the region of parameter 
space under consideration. 

5 Conclusions 

We have presented a novel solution to the dark matter (DM) problem, provided by a weakly 
interacting composite massive particle. This solution is based on using composite pseudo 
Nambu-Coldstone bosons not to solve the hierarchy problem but to provide a good DM 
candidate. Indeed, as depicted in figure 5, once we shift the focus from the hierarchy 
problem to the DM explanation (allowing the Higgs boson to be fully elementary too), 
then the minuses in Composite Higgs Models (CHMs) become pros in our new setup. In 
particular, CHMs require the Standard Model (SM) fermions to mix linearly with fermionic 
resonances in order to achieve the electroweak symmetry breaking (EWSB). It is well 
known that this mechanism can only provide a light Higgs boson in minimal models if 
the composite resonances are < TeV, which is more and more in tension with the current 
LHC data. As a matter of fact, the latest searches of vector-like quarks [22] exclude top 
partners with masses below 750 CeV in a complete model-independent way, and even 
masses around ~ 950 CeV can be reached depending of the decay mode. This is not longer 
true in Composite Dark Sectors, where the DM particles are allowed to be heavy enough 
to account for the relic density measurement, thus setting the composite resonances > few 
TeV and hence not in conflict with the current LHC data. Besides, the masses (and in 

^®Similar conceptual approaches have been previously considered in [10-13]. 
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Figure 5. Schematic representation of the Composite Higgs model versus the Composite Dark 
Sector formalism. 


fact the whole potential) of the scalar DM particles are only induced by loops of SM gauge 
bosons, what guarantees the DM stability even after EWSB. 

In order to quantitatively discuss all these features, we have considered in detail two 
minimal realizations of Composite Dark Sectors, corresponding to the cosets [SU{2)'^ x 
U{1)]/[SU{2) X U{1)] and SU{3)/[SU(2) x U{1)], which give rise to a real scalar triplet 
with y = 0 and a complex scalar doublet with Y = 1/2, respectively. In section 3 we have 
worked out the gauge interactions as well as the induced Coleman-Weinberg potential. 
For such a purpose, we consider a dual description of the composite sector with modified 
boundary conditions with respect to CHMs. We have explicitly shown that the minimum 
of the DM potential is always aligned in the direction that preserves the EW symmetry. We 
have also computed the masses of the neutral and charged scalars and the composite vector 
resonances as a function of the only two free parameters in the model, namely fn and go, 
corresponding to the typical scale and coupling of the composite sector. These two simple 
realizations turn out to correctly describe the DM phenomenology. We have considered 
experimental constraints from relic density measurements, direct detection experiments. 
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EW precision data —including the effects of both the scalars and the heavier resonances on 
the oblique parameters—, LHC searches of long-lived charged particles, collider searches of 
heavy narrow resonances —dijet, tt and dilepton final states— and LHC searches of monojet 
events. All together, they bound a large region of the fu — go plane, while allowing the 
dark sector to be the single component of the DM content in the universe (see figure 4). 
As it has been previously pointed out (see e.g. [79]), there is a significant complementarity 
between collider and non-collider searches, that will certainly allow us to probe the complete 
parameter space region in near future experiments. In addition, these models can be 
disentangled from other DM solutions. In particular, two main predictions of our scenario 
are the presence of long-lived charged scalars in the mass region 300 GeV < m^± < 2000 
GeV and heavy narrow vector resonances decaying equally into all the SM fermions. These 
signatures contrast with those of other related models of DM previously considered in the 
literatute [ 6 , 8 , 12 , 13, 80-85]. A dedicated study of the reach of future experiments to 
distinguish among these different alternatives is however beyond of the scope of this paper, 
and will be considered elsewhere. 
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A Gronp theory of [SU{2)‘^ x U{1)]/[SU{2) x U{1)] 


Since the pNGBs are not charged under the unbroken 1/(1) (identified with the hyper¬ 
charge), the relevant group theory is actually encoded in SU{2) x SU{2). The three gen¬ 
erators of the left SU ( 2 ) can be chosen to be 
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/ 0 0 i 0 \ 
0 0 0 -i 
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while the right SU(2) is generated by 
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Indeed, for any i,j, [5');^,S']j] = 0 and where eijk stands for the 

totally antisymmetric tensor. The unbroken subgroup SU{2)y is generated by the linear 
combinations T* = 5^ -|- The electromagnetic operator is then given by Q = 5 * 3 . The 
coset generators are written as A* = 
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B Group theory of SU{3)/[SU(2) x f/(l)] 

In standard notation, the SU (3) generators are written as 
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^0 1 0 ; 



The matrices = S'^ and generate the unbroken SU{2) subgroup, while 

t/(l) is generated by S'®/\/3. The electromagnetic charge is thus written as Q = S®+S®/\/3. 
Finally, the coset generators are given by = S^,X^ = S^,X^ = S® and = S’^. 
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